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By varying the density of an ultracold 88 Sr sample from 10 9 cm -3 to > 10 12 cm -3 , we make 
the first definitive measurement of the density-related frequency shift and linewidth broadening of 
the 1 5o - 3 Pi optical clock transition in an alkaline earth system. In addition, we report the most 
accurate measurement to date of the 8S Sr 1 So — 3 Pi optical clock transition frequency. Including a 
detailed analysis of systematic errors, the frequency is (434 829 121 312 334 it 20 sta t ± 33 3!/s ) Hz. 

PACS numbers: 42.62.Eh; 32.80-t; 32.70.Jz; 42.62.Fi 
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Neutral-atom-based optical frequency standards are 
becoming serious contenders for the realization of highly 
stable and accurate optical atomic clocks 0, Li L| 
Progress in this field is rapid, with significant ad- 
vances in experimental techniques from several related 
fields including laser cooling to ultracold temperatures 
QLULEIIlIi trapping configurations suitable for highly ac- 
curate clocks Q , the development of high quality optical 
local oscillators 8], and precise frequency measurement 
and distribution enabled by femtosecond optical combs 

Several systems have been proposed for neutral-atom- 
based optical clocks, each emphasizing narrow intercom- 
bination transitions in alkaline earth pi I2I fl^L fill fl4| and 
Yb atoms [lH . High resolution spectroscopic tech- 
niques employed to date include recoil-free absorption 
inside 1-D optical lattices 0], free space Ramsey inter- 
rogation and atom interferometry |l7j , and free space sat- 
urated absorption [Tsj . The fermionic strontium isotope, 
87 Sr, offers a weakly allowed 1 So - 3 Po optical clock tran- 
sition with a potential resonance quality factor exceeding 
10 15 and reduced polarization dependence and collisional 
shifts of the clock frequency when the atoms are confined 
inside magic- wavelength optical lattices [3J, |20j . How- 
ever, the large nuclear spin (I = 9/2) brings complex- 
ity in state preparation and field control. Recently new 
schemes have emerged that take advantage of the most 
abundant spin-zero isotope ( 88 Sr) by engineering a 'So - 
3 Po clock transition with a perfect scalar nature pi) . 

In this Letter, we present precision spectroscopy of 
ultracold 88 Sr in free space. With reference to the Cs 
primary standard and with a detailed investigation of 
systematic frequency shifts, we determine the absolute 
frequency of the 1 Sq — 3 Pi clock transition with a sta- 
tistical (systematic) uncertainty of 20(33) Hz. This level 
of measurement precision represents an improvement of 
more than 200 times over recent studies performed with 
thermal atomic beams [22|. Collision-related frequency 
shifts can significantly impact both the stability and ac- 
curacy of microwave and optical frequency standards. 
These effects have been studied in detail for both Cs and 
Rb microwave clocks [23, 0, E3 • For optical frequency 
standards based on intercombination transitions in al- 
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FIG. 1: (color online) Top view of the ultracold Sr spec- 
trometer. The two counter-propagating beams are collimated 
and overlapped. The beam intensities are balanced via atomic 
fluorescence signals. Light polarization is parallel to the bias 
magnetic field. A simplified 88 Sr level diagram and relevant 
experimental parameters are given in the bottom two panels. 



kaline earth atoms, it is expected that the dependence 
of the fractional frequency shift on the atomic density 
is reduced by more than three orders of magnitude in 
comparison to Cs clocks j2(|. Previous efforts to mea- 
sure this effect in Ca, however, were hampered by low 
sample densities ~ 10 9 cm -3 , resulting in measured shift 
coefficients with errors exceeding their associated abso- 
lute values H3- To overcome this limitation, we have 
performed absolute frequency measurements with ~1 fiK 
88 Sr atoms for densities spanning three orders of magni- 
tude, from 10 9 cm~ 3 to 10 12 cm~ 3 . For the first time 
we report definitive density-related line center shifts and 
spectral broadening of an optical clock transition in ul- 
tracold atoms. We find the density-dependent fractional 
frequency shift is 250 times smaller than that of Cs. 

Figure 1 shows a simplified view of the ultracold 88 Sr 
spectrometer. Components pertainingto atom cooling 
and trapping are described elsewhere |7l |28|. Measure- 
ments reported here utilize 10 6 to 10 8 atoms cooled to 
~ 1.3 /iK. Spectroscopic probing of the 1 Sq - 3 Pi clock 
transition is performed in free space after the trapping 
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laser beams and the quadrupole magnetic field are ex 
tinguished. The atomic density n during spectroscopic 
probing is varied by changing either the initial atom num 
ber in the trap or the free-flight time after the atom 
are released to free space. The typical free-flight time, 
range from 5 to 20 ms. Spatial densities are determinec 
from calibrated fluorescence images. The probe lase 
linewidth is <100 Hz. The laser's absolute frequency i 
measured continuously via a femtosecond optical coml 
that is locked to a Cs-referenced hydrogen maser via i 
fiber optic link to NIST [3l]]. The probe laser is spli 
into two counter-propagating beams, each coupled re 
spectively into a single-mode, polarization-maintainini 
optical fiber before delivery to the atomic cloud. Tb 
probe beams have 1/e diameters of 3 mm, > 50 m ra 
dius of curvature at the atom cloud, and are precisel; 
overlapped (<50 zxrad mutual tilt angle) by back cou 
pling each beam into the opposing beam's fiber launcher 
Both beams are normal to gravity to within 1 mrad. T< 
selectively probe the 1 S'o (mj = 0) to 3 Pi (m'j, = 
transition, the linear probe beam polarization is paralle 
(within 10 mrad) to a 4.6(2) xlO -4 T bias magnetic fiek 
turned on during the measurement cycle. The counter 
propagating beam intensities are balanced to within 2.59 
by monitoring atomic fluorescence signals. The sing! 
beam intensity (steady-state saturation of ~ 1.1) corre 
sponds to a resonant 7r-pulse time of ~ 87 ixs comparec 
to the probe pulse duration of 100 /is. Each point in i 
given trace of the atomic resonance is independently nor 
malized by the simultaneously measured atom number. 

Before proceeding further it is important to emphasize 
that the counter-propagating beam configuration is es 
sential to obtaining accurate values for the 1 S'o - 3 Pi lin 
center frequency. The ultracold 88 Sr temperature offers 
an advantage in reducing systematic frequency shifts as- 
sociated with atomic motion. However, we find that due 
to trapping-beam-intensity imbalance the atomic cloud 
acquires typical drift velocities along the probe beam axis 
of ~1 mm/s after release from the trap |29(. Although 
small, this drift corresponds to a frequency shift of 1.5 
kHz if a single probe beam is used for the spectroscopy. 
This frequency offset is canceled without perturbing the 
line shape if the atomic sample is illuminated by two 
equal intensity, counter-propagating probe beams. The 
presence of both beams also enables cancellation of the 
first-order gravity- induced Doppler shift. For the more 
realistic ± 2.5% intensity balance achieved here, the line 
shape is unperturbed but the true line center remains 
uncertain by ± 20 Hz, an effect that could be eliminated 
in future measurements by using an optical cavity to es- 
tablish the two probe beams. Although the two beam 
configuration described here resembles traditional satu- 
rated absorption, the saturation effect is greatly reduced 
since the excitation pulse length is comparable to the 
spontaneous decay time (21 fjs) while, for the ultra-cold 
temperatures employed here, the homogeneous and in- 
homogeneous linewidths are comparable. We estimate 
the saturation effect causes a slight reduction of ~ 1% of 
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FIG. 2: (color online) Dependence of the 1 So — 3 -Pi transition 
linewidth on n for 10 9 cm -3 < n < 10 12 cm . (a) Resonance 
profiles obtained for two different n. Circles represent data 
and the solid lines are fits to Voigt profiles. f c , wl, and wg 
are the line center frequency, the Lorentzian linewidth, and 
the Gaussian linewidth, respectively, (b) Wl and wg vs n, 
indicating a constant wg and a linearly broadened wl- The 
inset is a semi-log plot of the same data. 



the fluorescence signal peak at (u>q + ur), where u>r is 
the one photon-recoil frequency (27r 4775 Hz) and ujq the 
original atomic resonance. At the location of stimulated 
emission at (uiq - ojr), the reduction of fluorescence sig- 
nal is estimated to be < 0.1%. Such a small deviation 
from the originally symmetrical lineshape causes a sys- 
tematic shift in the measured resonance frequency by ~ 
+4(20) Hz. Note also that Ramsey interrogation is not 
well suited to free-space measurements of the 88 Sr 1 S$ 
- 3 P\ transition as its linewidth is slightly larger than 
the single-photon recoil frequency. Hence the Ramsey 
time required to align the recoil components exceeds the 
excited state lifetime. 

Figure 2 illustrates the dependence of the 1 Sq - 3 Pi 
transition linewidth on the peak density n for n ranging 
from 10 9 cm -3 to 10 12 cm -3 . Figure 2(a) shows two typ- 
ical resonance profiles obtained at densities of 8.4 x 10 9 
cm~ 3 and 6.3 x 10 11 cm~ 3 . The frequency axis is cali- 
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FIG. 3: Dependence of the 1 S'o — 3 Pi transition line center on 
n. Filled circles (the solid line) are data (a linear fit). A frac- 
tional frequency shift of —2.9(7) x 10~ 24 cm 3 is determined. 
The inset is a semi-log plot of the same data. 



brated with respect to the Cs standard. Circles represent 
experimental data while the solid lines are fits to Voigt 
profiles. The fits establish the line-center f c , the Gaus- 
sian linewidth wq, and the Lorentzian linewidth wl- As 
n increases by a factor of ~ 80, both f c and wl change 
significantly beyond the associated measurement uncer- 
tainties. As expected, however, wq remains essentially 
constant across the entire density range, as confirmed 
by the flat linear fit to the open circles in Fig. 2(b). 
The measured wl, shown as filled circles in Fig. 2(b), 
demonstrates a significant linear broadening coefficient 
of 2.8(3) x 1CP 8 Hz cm 3 . The zero-density extrapolated 
linewidth of 14.5(5) kHz (Fig. 2(b) inset) agrees with 
the predicted linewidth of 13.6(5) kHz after consider- 
ing power and interaction-time broadenings. In addition, 
we find that the integrated area under the Voigt profile, 
which is proportional to the number of photons per atom 
that contribute to the observed fluorescence signal, re- 
mains constant within the measurement uncertainty for 
n < 5 x 10 11 cm~ 3 . For n > 5 x 10 11 cm -3 (peak optical 
density of ~ 4.5), however, the Voigt profile area de- 
creases steadily with n, an effect that arises from probe 
beam attenuation and possibly radiation trapping. 

Figure 3 summarizes the frequency shift of the So - 
3 Pi transition vs n for n < 1 x 10 12 cm -3 . Filled cir- 
cles (the solid line) are experimental data (a linear fit). 
From the fit we determine a fractional frequency shift of 
-2.9(7) x 1(T 24 cm 3 , ~250 times smaller than the den- 
sity shift for Cs [24[. For n > 10 12 cm~ 3 (not shown), 
observed frequency shifts deviate significantly from the 
linear dependence displayed for n < 5 x 10 11 cm -3 . Re- 
solving this deviation requires further investigation. 

The density-related frequency shift and linewidth- 
broadening coefficients described above both exceed pre- 
dicted values based on general S-matrix calculations of s- 
wave collisions [2(| [3(j ■ According to Ref . [2(| , the clock 
shift (in Hz) is — An(h/m)L s , where m is the atomic 
mass, 2ttH is Planck's constant, and L s is a length that 
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FIG. 4: Accumulated record of the zero density 1 S'o — 3 
Pi optical clock transition frequency. Each entry is a 
weighted mean of several independent measurements that 
are referenced to the NIST Cs fountain clock. The over- 
all weighted mean value is (before correction for systematic 
shifts) 434 829 121 317 113 Hz, with a statistical uncertainty 
of 20 Hz. 



for finite collision temperatures is bound to the order of 
l/k. Here hk = \pmE is the thermal collision momentum 
and E is the collision energy. The bound for the magni- 
tude of the S-matrix is unity. The maximum fractional 
frequency shift is thus predicted to be —0.5 x 10 -24 cm 3 , 
^l/6 th of the measured shift. In the low temperature 
limit and in the absence of resonant scattering, L s is di- 
rectly related to the scattering lengths of the two atomic 
states associated with the clock transition. 

Elastic collisions do not account for the relatively large 
linewidth (wl) broadening observed here. For inelastic 
contributions, quenching of the 3 Pi excited state was ex- 
pected to be dominated by 3 Pi + 3 P\ rather than 3 Pi + 
1 Sq collisions. The quenching rate constant for 3 Pi + 3 Pi 
collisions has not been determined previously. An esti- 
mate based on the S-matrix formalism described above, 
however, limits the linewidth broadening coefficient to 
~2x 10~ 10 Hz cm 3 , or 140 times smaller than the mea- 
sured value. Other mechanisms are then most likely re- 
sponsible for the observed linewidth broadening. This 
conclusion is supported by the experimental observation 
that changing the fractional 3 Pi density does not signifi- 
cantly alter wl, indicating 3 Pi + 3 Pi quenching collisions 
are not the primary broadening mechanism. We have 
also experimentally ruled out any significant linewidth 
contributions due to depolarization of 3 Pi to the lower 
Zeeman level by 1 Sq collisions in the presence of the bias 
magnetic field. Note, however, that the relatively large 
Sr C3 coefficient leads to non-negligible p-wave contribu- 
tions to the 3 Pi + 1 Sq process even at a temperature 
^400 nK. Thus higher order partial waves may make 
significant contributions to the scattering process. An- 
other experimental observation reveals that for a given 
range of n, a longer free-flight time (compensated by us- 
ing a larger initial number of atoms) gives rise to a larger 
broadening coefficient. Here it is likely that the effec- 
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TABLE I: Systematic corrections and their associated uncer- 
tainties for the absolute frequency of the 1 S'o - 3 -Pi clock tran- 
sition. 



Contributor 


Correction (Hz) 


Uncertainty (Hz 


mutual beam tilt 





5 


position offset 





2 


velocity 





4 


2nd order Doppler 





2xl0~ 4 


1st order Zeeman 





0.02 


2nd order Zeeman 


4 


0.3 


light shift 





0.01 


blackbody shift 


1 


0.1 


recoil shift 


-4775 


< lxlO" 3 


stimulated emission dip 


-4 


20 


probe power balance (I) 


-5 


15 


probe power balance (II) 





20 


density shift 





2.9 


maser calibration 





1.7 


Systematic totals 


-4779 


33 



tive thermal collision momentum k among neighboring 
atoms is reduced due to their directional correlations in 
the free-flight expansion, resulting in a larger bound for 
L s since the bound is set by 1/k. We note that the linear 
fit in Fig. 2 includes the entire data set taken at various 
free- flight times (5 ms the shortest time). 

Figure 4 shows a record of the absolute 1 S'o — 3 Pi opti- 
cal frequency accumulated over two months. Each entry 
represents the weighted average of multiple independent 
measurements. A majority of the measurements are per- 
formed in the low density regime (n < 10 10 cm -3 ) and all 
measurements are corrected for the density shift shown in 
Fig. 3. The final weighted mean is 434 829 121 317 113 
Hz, with a statistical uncertainty of 20 Hz. A detailed 
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list of the systematic frequency shifts and their uncer- 
tainties are given in Table I. The first four items detail 
the residual Doppler effect from the atomic motion in 
free space. The next four items describe frequency shifts 
due to residual electric and magnetic fields. All have 
negligible contributions to the present uncertainty. The 
photon recoil correction is known to 1 mHz. To ensure 
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center occurs for a factor of two decrease in the overall 
optical power. Lineshape-asymmetry-related frequency 
shift caused by the stimulated-emission photon recoil is 
^4(20) Hz. The other dominant systematic uncertainty 
arises from the probe beam intensity imbalance. Along 
with the 20 Hz uncertainty described previously (Table 
1, probe power balance II), intensity mismatch leads to 
preferential multiple scattering in favor of the stronger 
beam, resulting in a blue shift due to photon recoil of 
5(15) Hz (probe power balance I). The overall correction 
to the measured frequency is therefore —4779 Hz, with a 
conservatively estimated systematic uncertainty of ~ 33 
Hz. Applying this correction, the 1 Sq — 3 P\ transition 
frequency is (434 829 121 312 334 ± 20 stat ± ii svs ) Hz. 

In summary, we have presented the most accurate 
frequency measurement of a Sr-based optical frequency 
standard with a detailed investigation of possible system- 
atic sources of error. Using the dense, ultracold sample, 
we have also performed the first definitive measurement 
of density-related frequency shifts and linewidth broad- 
ening for an optical clock transition. 
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